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The Stokes V parameter characterizes asymmetry of amplitudes between right- and left-handed 
waves, and non-vanishing value of the V parameter yields a circularly polarized signal. Cosmolog- 
ically, V parameter may be a direct probe for parity violation in the universe. In this paper, we 
theoretically investigate a measurement of this parameter, particularly focusing on the gravitational- 
wave backgrounds observed via ground-based interferometers. In contrast to the traditional anal- 
ysis that only considers the total amplitude (or equivalently J7gw), the signal analysis including a 
circular-polarized mode has a rich structure due to the multi-dimensionality of target parameters. 
We show that, by using the network of next-generation detectors, separation between polarized and 
unpolarized modes can be performed with small statistical loss induced by their correlation. 

PACS numbers: 



I. INTRODUCTION 

Because of the extremely weak signal, a direct detection of gravitational waves is a technically challenging issue, and 
we have not yet succeeded the direct detection despite extensive efforts. Nevertheless, the weakness of the gravitational 
interaction may be a great advantage for astronomy and cosmology, because gravitational waves can propagate to 
us from very early universe almost without scattering and absorption 0,0. In this respect, stochastic background 
of gravitational waves is one of the most important targets for gravitational wave astronomy [3|]. If detected, the 
stochastic background will serve as an invaluable fossil to study the physics at extremely high-energy scale for which 
other methods cannot be attainable. 

Over the last decade, sensitivity of gravitational wave detectors to the stochastic background has drastically im- 
proved. We will soon reach at the sensitivity level f^GW % 10~ 5 around 100Hz where Ogw is the energy density 
of the gravitational waves normalized by critical density of the universe. This level is below the indirect cosmological 
constraints, such as derived from the observed abundance of light elements 1 (see @ for the constraints from cosmic 
microwave background), and in this sense, gravitational wave detectors will provide a unique opportunity to directly 
constrain the early universe. In order to further get a stringent constraint and/or valuable information from the next- 
generation detectors, one important approach is to improve the statistical analysis of gravitational wave backgrounds. 
So far, most of theoretical studies on the gravitational-wave backgrounds have been directed to its energy spectrum 
(for its anisotropies, see, e.g., 0). The authors recently provided a brief sketch for measurement of the Stokes V pa- 
rameter of the gravitational-wave background via correlation analysis of ground-based detectors 0] (see [|[ for cosmic 
microwave background and [j| for space gravitational wave detectors such as LISA [T(3] or BBO/DECIGO (TTI. IT3|). 
The Stokes V parameter may be basic observable to quantify the parity violation process. One of such parity violation 
process is through the Chern-Simons term that might be originated from string theory [l3| . This paper is a follow-up 
study to the preceding short report. In addition to detailed explanations and supplementary materials to the previous 
paper, we developed a new statistical framework to deal with multiple parameters of gravitational wave backgrounds, 
and we specifically applied it to simultaneous estimation of amplitudes of both the unpolarized and polarized modes 
of the gravitational-wave background. 

This paper is organized as follows. In section [III we describe polarization decomposition of a gravitational-wave 
background, and define its Stokes V parameter. The basic framework to treat polarized gravitational waves is 
essentially the same one as in the case of electromagnetic waves [l5| . In section IIII1 we explain the correlation 
analysis for the gravitational-wave background and introduce the overlap functions that characterize sensitivities 
to the polarized and unpolarized modes. Then, we discuss basic properties of the overlap functions, and calculate 
them for the planed next-generation detectors, such as advanced LIGO. In section IIV1 broadband analysis of the 
gravitational-wave background is studied, taking into account the measurement noises for each detector. In section 
IVT we discuss how well we can separately measure the polarized and unpolarized modes. In contrast to the traditional 
arguments only for the unpolarized mode, the situation considered here is more complicated. We provide a statistical 
framework to analyze multiple parameters of the stochastic background with correlation analysis. Finally, section IVT1 
is a brief summary of this paper. Appendix [X] presents the derivation for the expressions of the overlap functions. 
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This geometrical derivation is similar to that given in Ref . [l8j . In appendix [B] we discuss the probability distribution 
functions (PDFs) of basic observational quantities with correlation analysis. In appendix [UJ we derive the formal 
expressions for optimal signal-to-noise ratios for detectors more than four. In appendix[Dl we comment on the surface 
of the Moon as potential sites for laser interferometers. 



II. CIRCULAR POLARIZATION 



Let us first describe the polarization states of stochastic gravitational waves. We consider a plane wave expansion 
of gravitational-wave backgrounds as 

h ij (t,x)= Yl r df I dn hp(f,n)e 2 ^- t+n - x ^(n). (1) 

P=+ x ^ 2 

Here, the bases for transverse-traceless tensor e p (P = +, x) are given as 

e + = eg ® eg ~ e^ ® &</,, e x = eg ® + ® eg (2) 

with the unit vectors eg and being normal to the propagation direction n that are associated with a right-handed 
Cartesian coordinate: 

eg = (cos cos</>, cos 8 sin 0, — sin 0), = (— sin 0, cos 0, 0). (3) 

On the other hand, the random amplitude hp represents the mode coefficients and the statistical properties of it are 
characterized by the power spectral density given by (hpi(n)h P2 (n')) (PI, P2 = +, x) for two polarization modes as 

'h + (f,n)h* + (f,n')} (h + (f,n)h* x (f',n'))\_l 2 l , ( I(f,n) + Q(f,n) U{f,n)-iV(f,n)\ 

[h x (f,n)hl(f,n>)) (h x (f,n)hUr,n>))) ~ 2° B[n " j ° DU 1 > \ U(f, n) + iV(f, n) I(f,n)-Q(f,n))> 

with delta functions <5d( - ) and the notation (• • •) for an ensemble average. Here, the quantities /, Q, U and V are the 
Stokes parameters and are real functions of direction n. Alternative to the linear polarization bases (e + , e x ), we may 
use the circular polarization bases (e R ,e L ) (right- and left-handed modes) 

r _ (e+ +ie x ) t _ (e+ -ie x ) 

for the plane wave expansion ((T|). Two coefficients hp.L for the corresponding modes are given as 

{h + -ih x ) (h++ih x ) 
kR= V2 ' kL= V2 ■ (6) 

Then the covariance matrix is recast as 

(h R (f,n)h R (f',n')*) (h L (f,n)h R (f,n'r)\_K ( „ „,y x ( f f n f I(f,n) + V(f,n) Q(f,n)-iU(f,n)\ 
(h R (f,n)h L (f,n>y) (h L (f,n)h L (f,n>y)) ~ 2 dD(n ^d(/ / ) \^ Q ^ n ) +i U(J,n) I(f,n)-V(f,n) J 

(7) 

With this expression, it is apparent that the real parameter V characterizes the asymmetry of amplitudes between 
right- and left-handed waves, while the parameter /(> |V|) represents their total amplitude. For example, if we can 
observationally establish V > 0, the background is dominated by right-handed waves. Since the parity transformation 
interchanges the two polarization modes, the asymmetry is closely related to parity violation process (see e.g. [l3l[l4j 
for recent theoretical studies). Therefore, we may detect signature of parity violation in the early universe by analyzing 
the V parameter of gravitational-wave backgrounds. This is the basic motivation of this paper. 

Since the two parameters / and V have spin 0, their angular dependence can be expanded by the standard (scalar) 
spherical harmonics Yi m : 

oo t oo i 

I(f,n)=J2 ^m(f)Y em (n), V(f,n) = Y^ E V tm (f)Y tm (n). (8) 

1=0 m=-i £=0 m=-i 
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On the other hand, the combinations Q±iU describe the linear polarization and have spin ±4 reflecting spin-2 nature 
of gravitational waves. They are expanded with the spin- weighted spherical harmonics as 

oo £ oo £ 

(Q + i U)(f,n) = E E PL(fUYUn), (Q - i U)(f,n) = E E P fm(fU Y tm(n). (9) 

£=4m=-£ £=4m=-£ 

Note that Q ±iU do not have monopole components (I = 0), because the linear modes introduce specific spatial 
directions. Since the observed universe is highly homogeneous and isotropic on large spatial scales, it is reasonable 
to set the monopole modes of a cosmological stochastic background as our primary targets. Therefore, in this paper, 
we do not study the linear polarization Q ±iU. From the same reason, we also neglect the directional dependence of 
the / and V modes. 

Next, we describe the frequency dependence of the gravitational-wave background. To characterize the gravitational 
waves in the cosmological context, rather than the spectral density, the normalized logarithmic energy density of the 
stochastic background, f2Gw(/), is frequently used in the literature [lalla]- The density Ogw(/) is defined by the 
spectral density / as 

4tt 2 f 3 

"gw / = —^-I(f), (10) 

Pc 

where p c (= 3Hq/8tt, Ho = 70/i7okm/sec/Mpc: the Hubble parameter) is the critical density of the universe. We also 
define the polarization degree by n(/) = V(f)/I(f). In terms of the quantities Hgw(/) and n(/), the asymmetry 
parameter V is expressed as 

= ^3^w(/)n(/). (li) 



III. OVERLAP FUNCTIONS FOR GROUND-BASED DETECTORS 



This section discusses the overlap functions for correlation signals as the basic ingredient for correlation analysis of 
gravitational-wave background. In Sec lIII Al the definition and the analytic formula for overlap functions are given. 
Subsequently, Sec lIII B | MI CI and IIII Dl discuss special or limiting cases for overlap functions in order to understand 
their geometrical properties. After describing some mathematical properties in Sec lIII E"l we evaluate the overlaps 
functions for specific pairs of five detectors in Sec lIII Fl 



A. Formulation 



Let us begin by considering how we can detect the monopole components of the I and V modes with laser interfer- 
ometers. Response H a of a detector a at x a is written as 

H a (f)= f dn V h P (f,n)F^(nJ)e 2 ^ n - x -. (12) 

Js2 P=+,x 

The function is the beam pattern function and it represents the response of the detector to each linear polarization 
mode. Here we used the conventional linear polarization bases. 

To distinguish the background signals from detector noises and to obtain a large signal-to-noise ratio, the correlation 
analysis with multiple detectors is essential [lg, [l7|, [l8|, [l9( . We define the correlation C a b(f) of data streams obtained 
from two detectors a and b as 

(H a (f)H b (f)*) ee C ab (f)S D (f - /')■ (13) 
Keeping the monopole contribution only, its expectation value is written as 

CM = y h Iiab (f)I(f) + lV.a b {f)V(f)\ , (14) 

where 7/ is the overlap function and given by [HI, [l!| 

U, ab (f) = ^ J s2 dn [{F+F+* + F* F b x * } e^ nm ] , (15) 
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with rewriting x a — x b = Dm (D: distance, m: unit vector) and y = 2nfD/c. The variable y represents the phase 
difference at two cites a and b for waves with a propagation direction m. Similarly, the function "fv,ab{f) is given by 

lv. ab {f) = f / dn [i {F+F** - F* F b + * } e iynrn ] . (16) 

»7T J S 2 

Two functions 7/ and 7y are purely determined by relative configuration of two detectors. 

Here, we summarize the response of ground-based L-shaped interferometer a. We assume that two arms of the 
next-generation interferometer have equal length with opening angle of 90°. We denote the unit vectors for the 
directions of its two arms as u and v. At the frequency regime where the wavelength of the incident gravitational 
wave is much longer than the armlength, the beam pattern function takes a simple form as 

Ff = d a : e p (n), (17) 

where the colon represents a double contraction and the detector tensor d a is given by 

(u a <g> U a - V b (gi v b ) 
da - ^ ' ^ ) 

In reality, there might be some exceptional cases that the opening angle of two arms is slightly different from 90° such 
as GEO600, whose opening angle is 94. 3° 12011 . However, the response of such a detector can be treated as the one of 
the right-angled interferometer (see eg. [HI for the case with LISA). 
Note that from equations (JTTJ) and (fT5|) , we have 



and Schwartz inequality implies 



s I ^W^-s /,*•['"''"] -? < 1<J) 

- 1 < 7i,oi < 1, - 1 < lv,ab < 1. (20) 

In Table|lJ we list positions and orientations of the ongoing (and planned) kilometer-size interferometers, AIGO 21], 
LCGT [22|], LIGO-Hanford, LIGO-Livineston J23I . and Virgo [24| . As a reference, we also list two sub- kilometer-size 
interferometers, TAMA300 [25| and GEO600 [20]. We use a spherical coordinate system {6, 4>) with which the north 
pole is at 8 = 0°, and <j> represents longitude. The orientation a is the angle between the local east direction and 
the bisecting line of two arms of each detector measured counterclockwise. Since the beam pattern functions have 
spin-2 character with respect to the rotation of detector, meaningful information here is the angle a module 90°. In 
what follows, we mainly focus on the first five detectors in Table U with their abbreviations (A,C,H,L,V) and with 
R = i?E = 6400km for the radius of the Earth 1 , but in section IIV CI we also discuss the detectors placed on the 
Moon as an exceptional but interesting case (see also appendix |D| . 





e 





a 


AIGO (A) 


121.4 


115.7 


-45.0 


LCGT (C) 


53.6 


137.3 


70.0 


LIGO Hanford (H) 


43.5 


-119.4 


171.8 


LIGO Livingston (L) 


59.4 


-90.8 


243.0 


Virgo (V) 


46.4 


10.5 


116.5 


TAMA300 


54.3 


139.5 


225.0 


GEO600 


47.7 


9.8 


68.8 



TABLE I: Positions (9, (f>) and orientation angles a of detectors (in units of degree) on the Earth. 



For the monopole modes of the stochastic background, only the relative configuration of two detectors is relevant 
with the correlation C a b and we do not need to deal with their overall rotation. Therefore, without loss of generality, 
their configuration is characterized by the three angular parameters (/?, <7i, 02), shown in Figure [1] [l8[. Here, (3 is the 
separation angle between two detectors measured from the center of the Earth. The angle o\ (01) is the orientation of 



1 We use the roman V for Virgo detector and the italic V for the polarization mode. 
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a 




FIG. 1: Detector planes are tangential to a sphere. Two detectors a and 6 are separated by the angle j3 measured from the 
center of the sphere. The angles o\ and 02 describe the orientation of bisectors of interferometers in a counter-clockwise manner 
relative to the great circle joining two sites. 



the bisector of two arms of the detector a (b respectively) measured in counter-clockwise manner relative to the great 
circle connecting a and b. Their distance is given by D = 2i?E sin(/3/2) that determines a characteristic frequency 
/d = c/D for the overlap functions. Following [lj|, we define the angles 

A S (21) 
and the geometrical information for possible pairs made from the five detectors are summarized in Table [TXJ 





A 


C 


H 


L 


V 


AIGO (A) 


* 


70.8, 58.1, 31.4 


135.6, 53.7, 45.1 


157.3, 38.0, 2.08 


121.4, 19.2, 60.8 


LCGT (C) 


-0.61, -0.58, 0.81 


* 


72.4, 89.1, 25.6 


99.2, 42.4, 68.1 


86.6, 28.9, 5.6 


LIGO Hanford (H) 


-0.82, -1.00, -0.007 


1.0, -0.21, 0.98 


* 


27.2, 45.3, 62.2 


79.6, 61.8, 55.1 


LIGO Livingston (L) 


-0.88, 0.99, 0.15 


-0.98, 0.04, -1.0 


-1.00, -0.36, -0.93 


* 


76.8, 26.7, 83.1 


Virgo (V) 


0.23, -0.45, -0.89 


-0.43, 0.92, 0.38 


-0.43, -0.76, -0.65 


-0.29, 0.89, -0.46 


* 



TABLE II: Upper right: angle parameters (/3, 8, A) for each pair of detectors in units of degree. Lower left: numerical values 
(cos(45), cos(4A), sin(4A)) for each pair of detectors. 

The angular integral (|15[) can be performed analytically with explicit forms of the pattern functions, and we get 

li,ah = 6 1 (y,/3)cos(4 ( 5) + 9 2 (y, 0) cos(4A), (22) 

with 



<->!'.'/■ ('"]) (./■, • 1.12 • j^M 



(23) 



and 



3 45 169. \ A 5 27 \ „ / 1 5 3 , 

8 J0 + 56 J2 896 J V + {2 J °- 7 j2 ~ 224 j4 C ° S(3+ l~8 j0 " 56 j2 " 896 j4 ' (24) 



e 2 (y,/?) = 

The function j n is the n-th spherical Bessel function with its argument 



2nfD 47r/i? E . (8 

V = = sin - 

c c \ 2 



(25) 
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The expression ([22]) coincides with the formula (4.1) in Ref. [l8l |. 

In a similar manner, the overlap function for the V mode is given by 



with 



e 3 (y,/3) 



7v>& = e 3 (2/,/3) sin(4A) 
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sin 



-ji + sh ) + ( h + -33 ) cos (3 

and 



(26) 
(27) 

can be deduced from the 



Note that the dependence of the angles S and A on the overlap functions 
symmetric reasons Q. 

In appendix [XI we present a brief sketch to derive the expressions 7/ and 7y, using the symmetries of tensorial 
structure. Since our primary interest here is the dependence on the frequency / and the angle /?, we mainly use the 
set of the variables (/, (3), instead of (y, (3). 



B. Special cases and asymptotic profiles 

In order to get a physical insight into the overlap functions, it is instructive to consider geometrically simple 
configurations for two detectors. When a pair of detectors are placed on the same plane ((3 = 0°) and at the same 
position (D = 0), we have the identity (9i, 82) — (1, 0) and thus 7j )0 6 = cos(4<5). In contrast, for V mode, we obtain 
Iv.ab = for the coplanar configuration (/3 = 0°) and this is even true with finite separation D 7^ 0. The reason 
for this is explained in next subsection. Equation l|22p and the identity 82(2/, 0°) = indicates that the function 7/ 
depends very weakly on the parameter A at small angle (3 -C 180°. For ground-based detectors, the functions Qi(y, (3) 
depend on the angle (3 also through the variable y — 4ttRbJ ' sin(/3 /2) / c. Taking into account this fact, we obtain the 
following asymptotic profiles at small f3 (in unit of radian): 

e 1 = o(f3°), e 2 = o(/3 4 ), e 3 = o(/? 3 ). (28) 

On the other hand, for pair of detectors located at antipodal positions (/3 = 7r), the overlap function 7/ j0 b does not 
depend on the parameter S because of 0i(y, 180°) = 0. In this case, the asymptotic profiles become 

e 1 = o((^-/3) 4 ), e 2 = o((7r-/3)°), e 3 = o((7r-/3)°). (29) 

Note that at j3 = and tt, we have 

000! = 9^02 = 0003 = 0. (30) 



C. Coplanar configuration 

An L-shaped detector measures difference of spatial deformation towards its orthogonal two arms. This is purely 
geometrical measurement. If two detectors are placed on the same plane z = 0, there is an apparent geometrical 
symmetry for the system with respect to the plane. Due to the mirror symmetry to the plane, a right-handed wave 
coming from the direction (n Xl n yi n z ) and a left-handed wave from the direction [n x ,n y , —n z ), provide an identical 
correlation signal, if they have the same frequency and amplitude. Therefore, for an isotropic background, right- 
handed waves coming from two directions (n x ,n y , ±n z ) exactly cancel out in the correlation signal. The same is true 
for left-handed waves. As a result, the symmetric system has no sensitivity to the isotropic component of the V^-mode 
• We can directly confirm this cancellation from the definition (TT5)) and the following relations 

F+(n x ,n y ,n z ) = F+ (n x ,n y , -n z ), F* (n x ,n y ,n z ) = -F* (n x ,n y , -n z ), (31) 

which are easily derived from the symmetries of the polarization bases e + ' x (n) [28|. The cancellation of correlation 
signal is particularly important for setting orbits of space-based interferometers, such as BBO/DECIGO [9j. For 
detecting the monopole of the V^-mode, it is essential to break the symmetric configuration. 
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D. Optimal configuration 

In this subsection, we consider optimal configurations of two detectors (a, b) for measuring the I and V modes 
of stochastic backgrounds. To investigate the optimized parameters for overlap functions, there are two relevant 
issues; maximization of the signals ji^b an d Jv.ab, and switching off either of them (ji^b = or 7/ )0 & = 0) for their 
decomposition. For simplicity, we consider how to set the second detector b relative to the fixed first one a for a given 
separation angle 0. In this case, the sensitivities to the I- and y-modes are characterized by the remaining adjustable 
parameters, ct\ and 02- The former determines the position of the detector b, while the latter specifies its orientation 
(see FigfT]). Based on the expressions (|2"2")l and (j2"u) . one finds that there are three possibilities for the optimal detector 
orientation: 

Type I : cos(4A) = - cos(4<5) = ±1 (32) 

or 

Type II : cos(4A) = cos(4(5) = ±1 (33) 

to maximize the normalized SNR Sj^ a b [HI, and 

Type III : cos (4A) = cos (4(5) = (34) 

to erase the contribution from J- mode. The relative signs of the two functions Oi and 02 determine whether type I 
or type II is the optimal choice. 

For type I, the solutions of the two angles a\2 are a\ = a% = 45° (mod 90°) and the detector b must be placed 
on one of the two great circles passing through the detector a, parallel to one of the two arms as shown in Figure [21 
For a given separation /3, there are four points for the cites of the detector b. At each point we have four equivalent 
orientations as shown in the right panel of Figure [U This is because response of an L-shaped detector has mod-90° 
effective equivalence. After all, for a given separation j3, there are totally 4x4= 16 possible configurations of detector 
b. 

For type II, the second detector must reside in two great circles parallel or perpendicular to the bisecting line of 
each detector, as shown in Figure [31 As in the case of type I, with a given separation (3 we have totally 16 candidates 
for detector b. At the limits (3 — > 0° and (3 — > 180°, there are no essential differences between types I and II. 

Similarly, the type III configuration is realized by placing the second detector on one of the four great circles defined 
for types I and II, with rotating 45° relative to the first detector (see Fig. [3]). In this case we have 8 x 4 = 32 possible 
configurations for detector b. Note that the sensitivity to the V-mode is automatically switched off for the type I 
and II configurations and is conversely maximized for the type III configuration. This is because the function 7y is 
proportional to sin(4A). 

E. Basic properties of functions Oi 

In this subsection, specifically focusing on the detectors on the Earth with radius i?E=6400km, we study basic 
properties of the three functions 0i, 02 and ©3 in some details. Note that in general, for a sphere with radius i? s , 
there is one characteristic frequency c/R s and our results for the Earth at frequency / can be rescaled to those for 
the sphere with scaled frequency (Re/R s ) f 2 ■ Hence, the result presented here may be interpreted as the one for an 
arbitrary sphere, including multiple detectors placed on the Moon. 

In left panel of Figured! the function 0i(/, (3) is plotted against the angle parameter (3 at specific frequencies 
/ = 10, 50 and 70Hz. As shown in Sec. IIII CI we have 0i = 1 at ft = 0° that is the maximum value for 7/ for given 
frequency /. At frequency / > 10Hz relevant for ground-based detectors, the function |©i| becomes very small for 
a separation angle (3 > 90°, and we identically have 0i = at antipodal configuration (3 — 180°. In right panel of 
Figure SI the shape of the second function 02(/, (3) is shown. As discussed in Sec. IIII CI the function ©2 becomes 
vanishing at j3 = 0°. This function exhibits an oscillatory behavior in the range 0° < (3 < 180°, and the number of 
its nodes is approximately proportional to /-Re (see appendix [AJ . 

In Figure O we plot the overlap function |7/| for two optimal configurations, types I and II, at specific frequencies 
10 and 50Hz. The thin lines are for the type I with |7/| = |©i — ©2!, while the thick lines are for the type II with 



2 This is easily deduced from the fact that the functions Gi, 02 and O3 depend on frequency / only through the product /-Re 
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: b 
I 




Type I 

FIG. 2: Type I configuration with a given separation angle /3. Relative to a fixed L-shaped interferometer a, the second one 
must be placed on two great circles shown with long-dashed lines (left panel). We also have four equivalent detector orientations 
due to mod-90° freedom as shown in the right panel. 




Type I Type II Type III 



FIG. 3: Position and orientation of the second detector b relative to the fixed first one a. The long dashed lines are great 
circles passing the first one a. 



\ji\ = |@i + B2I. For angles close to (3 = 0° and 180°, two lines are almost identical. This is because only one 
component is dominant there, that is, |6i| 3> | ©2 1 at /? ~ 0°, and |0 2 | 3> |Oi| at (3 ~ 180°. Two components have 
comparable magnitude at j3 ~ 120° for / = 10Hz and at (3 ~ 60° for / = 50Hz. As shown with the curves for 
/ = 50Hz, both types I and II have chance to give the maximum value of |7/| for given (f,(3), depending on the 
relative signs of 61 and 62- 

Next, in Figure El the function 83 for the I^-mode is plotted. Note that we have |tv| = | ©3 1 for the type III 
configuration. As in the case for 82, the oscillating profiles have a number of nodes approximately proportional to 
/f?E- For given frequency /, we define the separation angle /3 ma x that maximizes the function | O3 1 in the range 
f3 e [0°, 180°]. In contrast to the simple results for the J-mode with max7/(/,/3) = 0i(/, (3 = 0°) = 1, the angle 
Anax defined for the V-mode is slightly complicated and it does depend on the frequency /. Figure [7] shows the angle 
(3max in unit of radian, plotted against the frequency. The frequency dependence in the range < / < 16.7Hz can be 
understood with the following three steps: 

(i) As commented earlier, we have dpQ^(f,(3 — 180°) = representing that the end point (3 = 180° is generally an 

extreme. At low frequency regime, the oscillating feature of 83 is relatively simple (see Fig. [HI, and the end 
point 8 3 (/,/3 = 180°) is the global maxima. We find /3 max = 180° for / < 12.8Hz. 

(ii) At / = 12.8Hz, the end point = 180° becomes an inflection point with <9|8 3 (/,/3 = 180°) = 0. Then, for 

/ > 12.8Hz, there appears a local maxima for 83 at f3 < 180° that determines the separation angle /3 ma x, as 
in Figure [JJ Meanwhile the end point 83(7, (3 = 180°) is now a local minimum. With increasing / it decreases 
and crosses at / = 15.7Hz. 
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FIG. 4: The functions Oi(/, (3) and 02(/, 13) for detectors on the Earth at frequencies / = 10Hz, 50Hz and 70Hz. 




FIG. 5: The optimal combinations |8i(/,/3) + 6i(/,/3)| (type II: thick lines) and |8i(/,/3) - Qi{f,f3)\ (type I: thin lines) for 
detectors on the Earth at frequencies / = 10Hz and 50Hz. 



(iii) The local maxima 6 3 (/,/? max ) at /3 max < 180° coincides with -6 3 (/,/3 = 180°) (> 0) at / = 16.7Hz, and the 
separation angle /3 max shows a discontinuous transition up to /3 max = 180° at / = 16.7Hz. 

We can observe similar cycles for the angle /3 max at / > 16.7Hz. The frequency dependent angle /3 max should 
be regarded as the optimal separation for narrow band detection for the y-mode signal. In Figure [5J we show the 
maximum value |9 3 (/, /? max )| as well as \@ 3 (f,0= 180°)| and |6 3 (/,/3 = 150°)|. The choice (3 = 150° is just for an 
example. The first two curves coincide at some frequency bands (as shown in Fig. [7] for /3 max ), while the example 
| @3 (/, (3 — 150°)| contacts with the dashed curve for maximum value |0 3 (/, (3 m &x\ only at specific discrete frequencies. 
In the two dimensional region with 0° < (3 < 180° and / > 0, the global maximum for 3 is 

- — (2cos2- 5sin2) = 0.84, (35) 
32 

which appear at (3 = 180° and / = c/2ttRe = 7.5Hz for detectors on the Earth with i?E = 6400km. In general, the 
function 3 is maximized at antipodal configuration {(3 = 180°) with / = c/2nR s , or equivalently y — 2 due to the 
scaling commented in the beginning of this subsection. Although, for practical purpose to detect the V-mode signal, 
the broad-band analysis with multiple detectors is essential, which we will discuss in next section, it is clear from 
Figure [8J that the separation (3 — 180° seems the best choice for detectors, whose bandwidths are much larger than 
the individual wiggle structure in this figure. 
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FIG. 6: The function 03(/, /3) for detectors on the Earth at frequencies / = 10Hz, 50Hz and 70Hz. 
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FIG. 7: Values of /3 max for detectors on the Earth as a function of frequency. While vertical lines are shown due to a software 
reason, there are discontinuities from /3 ma x < 180° to /3 max = 180°. 



F. Overlap functions for specific pairs of detectors 



Now, we analyze the geometry of ten pairs made from the five detectors listed in Table [T] In this paper, we do not 
consider the co-located and co-aligned pair of detectors, such as two interferometers (4km+2km) at LIGO-Hanford. 
Co-located and co-aligned detectors are possibly contaminated by the measurement noises which are statistically 
correlated with each other, making the detection of stochastic signals difficult. 

Let us first examine how well the pairs of existing or planned interferometers are suitable for I- and V-mode 
detection by comparing the angle parameters with those of the optimal configuration discussed in Sec. IIII Dl In left 
panel of Figure [H we plot the combination (cos(4<5), cos(4A)). ^,From this plot, the AL and AH pairs are found to 
be very close to the type I and type II configurations, respectively. Except for these, however, there are no other 
noticeable pairs. Turn to next consider a large separation angle (3 ~ 180°, where the parameter 8 becomes unimportant 
and the correlation signal can be approximately described by 

C ab y [63(2/, p) cos(4A) + Q 3 (y, 0) sin(4A)] . (36) 

Thus, in this case, the angle parameters (3 and A now play an important role. Since the regime (3 > 90° is preferable 
for the V-vaode detection, we next plot the combination (/3,cos(4A)) in right panel of Figure [51 Among various 
pairs of detectors, the CL pair realizes nearly ideal angle (sin(4A) = —1), although the separation angle of CL 
pair is intermediate, i.e., j3 = 99.2°. Other interesting pairs for the V mode with relatively large |sin(4A)| are AV 
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FIG. 8: The function -yv | for detectors on the Earth with type III configuration. The solid curve (dotted curve) is the result 
with f3 = 7r (P = 57r/6). The dashed line is result with I3 max for which the function |7v| becomes maximum with given frequency 
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FIG. 9: Left: Distribution of the combinations (cos 45, cos 4A) for ten detector pairs shown in Table ITT1 Points for three types 
of configurations are also given with squares. Right: Distribution of the combinations (/?, cos(4A)) for ten pairs shown in Table 
HTl At relatively large (3 the sensitivity to the V-mode is roughly proportional to sin(4A). 



(sin(4A) = -0.89) and CH (sin(4A) = 0.98). The HL pair has sin(4A) = -0.93, but its separation is small, /3 = 27.2°, 
where the amplitude 03(y,/3) is relatively small. 

Based on these considerations, in Figure [TUl we compile the overlap functions (77, 7v) for the ten pairs of detectors. 
There is characteristic frequency- width, A/ oc (sin/3/2) -1 , determined by the arrival-time difference of gravitational 
waves between two cites. The frequency interval is largest for the HL pair. For high frequencies, the peaks for the 
functions 7/ and jy have 1/4-cycle phase difference, as discussed in appendix [3] The AH pair is almost insensitive 
to the V mode, because it is close to the type II configuration. The situation is similar for the AL pair. Note that 
the CL and AV pairs have relatively good sensitivity to the V mode, as anticipated from the angular parameters. 
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FIG. 10: Overlap functions for specific pair of detector made from the five detectors, A, C, H, L and V. The solid lines are for 
7/ and the dotted lines for jv- 

IV. BROADBAND SIGNAL ANALYSIS 
A. Preliminary 

So far, we have only dealt with the correlation signal of gravitational-wave backgrounds. In practice, the signal 
is contaminated by detector's noises, and thus the broadband signal analysis is essential ingredient for detection of 
background signals with high signal-to-noise ratio. 

We model the data stream s a of a detector a by a summation of gravitational-wave signal H a and detector noise 
n a , 

s a = H a + n a . (37) 

Throughout this paper, we assume that the noise of detector, n a , obeys stationary and random processes and the 
noise correlation between any pair of detectors can be safely neglected. Then, covariance of the detector noises can 
be expressed as 

M/HCf')*) = ~<U£ D (/ - f'Wf), (38) 

where N a is the noise spectral density for detector a. 

To estimate the sensitivity of each pair of detectors, let us consider the simple case with correlation C a b(f) of two 
detectors a and b. As it has been shown in the literature, the total signal-to-noise ratio (SNR) is given by [H, [l9[ 
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(see also appendix |B|) 
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16vr 
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IOtt 2 
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(7/^ + 7v^) 2 



obs 



4f 



Kb(f) 

^gw(/) 2 (7/ 



- 1 



/W ab (/) 



(39) 



with the quantity Af a b defined by N a b = N a (f)Nf,(f). This is the result obtained in the weak-signal limit /(/) <C 
N{ a ,b\{f)- Note that the above formula just represents the SNR for the total amplitude of the background signals, 
and it does not imply the SNR for a pure I- or T^-mode signal. The separation of I- and T^-modes will be discussed 
in next section. 

For quantitative evaluation of SNRs, we need an explicit form of the noise spectral density. In the following, we use 
the fitting form of the noise spectra for advanced LIGO detector, N\i go . Assuming that all the detectors have identical 
noise spectra with Nn go , signal-to-noise ratios of stochastic signals are estimated. Based on Ref.[26[, the analytical fit 
of the noise spectrum -/Vij go is given by 



Mi go (/) 
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10Hz 

-46 / / 
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-47.25 
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-1.7 



1000Hz 
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for 240Hz < / < 3000Hz, 
otherwise. 



(40) 



The expression (|3"9")l implies that the weight function for SNR per logarithmic frequency interval d In / is proportional 
to (N(f)f 5 / 2 )^ 1 for flat input riGw(/)(7/+7vn) = const. In Figure ITTT using the analytic form of the noise spectrum, 
we plot the weight function. It becomes maximum around ~ 50Hz with its bandwidth ~ 100Hz. Note that the shape 
of the weight function for stochastic signals is close to the one for the signals produced by binary neutron stars, in 
which case the detectable distance, as the integral of weight function, is roughly proportional to 



dhxf 



1 



1/2 



(41) 



The next-generation detectors are primarily designed to have the good sensitivity to a chirping signal of binary 
neutron stars, and they are planned to achieve the similar performance for detecting these binaries. In this sense, our 
assumption that all the detectors have identical noise spectrum with advanced LIGO is reasonable. 
Finally, as a reference, we present the SNR for coincident detectors (77,06 = l,7v,ab = 0): 



SNR 



1/2 



^GW^70 



V 3yr J V 10-9 

This value will be frequently referred, as a baseline of the SNR for various situations considered below. 



(42) 



B. Signal-to- noise ratios for pair of detector 



The total SNR (|39|) depends strongly on model parameters of the background, including the polarization degree n. 
In order to present our numerical results concisely, we use the normalized form 

^ SNR 

^ SNRo' (43) 

To characterize sensitivity of each pair to I- or T^-mode signal, based on the above equation, we respectively define 
Qi and Qv by replacing {/, V} in the expression of SNR with p c /(AT: 2 f 3 ){D, GW , 0} and p c /(47r 2 / 3 ){0, Hgw}- 
These normalized SNRs Qi and Qv can be regarded as a rms value of overlap functions with the weight function 
(f 5/2 N(f))~K 

In Figure [T2l we present the normalized SNR for the optimal geometry, i.e., types I, II and III configurations 
(short-dashed, long-dashed and solid lines, respectively). One noticeable point is that a widely separated ((3 ~ 180°) 
pair is powerful to search for the V mode. At (3 — 180°, we find the following asymptotic relations 

5 . , . , sin y 5 .... , cos y . , , . 

7J ~ cos 4A y - , lv - - sin 4A y - 44 

2 y 2 y 
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FIG. 11: The weight function (N(f)f 5/2 ) for advanced LIGO. 

with y — AnRf jc. For detectors on the Earth, the characteristic frequency interval is c/2i?E ~ 20Hz, which is enough 
inside the bandwidth of advance LIGO, Af ~ 100Hz. As a result, the oscillation of the overlap functions are averaged 
out and the normalized SNRs Qj and Qv give a similar output for optimal configurations (types I-III) at f3 = 180°. 

In Figure [T3] and Table IIII1 we show the normalized SNRs for pairs made from the five interferometers in Table fl] 
To reduce the contribution from the J-mode signal, pairs that have been regarded as disadvantageous for constraining 
SJgw, can now play important roles for measuring the V mode, according to equation (|36p . The HL pair with 
cos(4(S) ~ 1 and sin(4A) ~ 0.93 realizes nearly maximum values simultaneously for Qi t ab arid Qv,ab, at its separation 
(3 = 27.2°. This is because Qi y ab is mainly determined by the angle 6 at a small (3, while Qv.ab depends only on A. 
This pair has the largest Qi among ten pairs of detectors. In contrast, the CL pair has good sensitivity to V, although 
it is relatively insensitive to the I mode, because of sin(4A) ~ 1 and | cos(4A)| = 0.04. Indeed, the orientation of the 
LCGT detector is only 1.2° different from the optimal direction (cos(4A) = 0) with respect to the LIGO-Livingston 
cite. As other interesting pairs, the AH pair is almost insensitive to the V mode with sin(4A) = —0.007 (nearly type 
II configuration with LIGO-Hanford) . The AV pair has a large Qv with | sin(4A)| = 0.89, but its Qi is much larger 
than that of the CL pair. 
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TABLE III: Normalized SNRs Qj (upper-right) and Qv (lower-left). 



C. Antipodal detectors 

So far, we have studied pairs of detectors on the Earth, strictly keeping the radius of sphere _Re = 6400km, for 
which the antipodal type III configuration is turned out to be optimal to realize the largest SNR Qv- Here, we discuss 
to what extent one can improve the sensitivity to the V mode by varying the radius of sphere. 

In Figure [HI we plot the broadband sensitivities Qi and Qv for antipodal detectors as functions of the radius 
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FIG. 12: Normalized signal to noise ratios (Qi, a b an< l Qv.ab) with optimal configurations for the I-mode (short dashed curve: 
type I, long dashed curve: type II) and for the V-mode (solid curve: type III with setting II = 1 for illustrative purpose). We 
use the noise curve for the advanced LIGO. 





FIG. 13: The normalized SNRs Qi (circles) and Qy (triangles) for detector pairs. 
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FIG. 14: The normalized SNR for antipodal configuration with radius R. Advanced LIGO noise curve is used. 
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FIG. 15: Overlap functions for antipodal detector pairs (solid curves: 71, dotted curves: jv)- 



of sphere, R. Note that the noise spectrum is fixed to Nu go (f) as before. Here, we put cos(4A) = 1 for Qj and 
sin(4A) = 1 for Qv- While the value Qi is maximized at R = and we obtain Qi — 1 and Qv — 0, the maximum 
value of Qv is achieved when R = 600km, leading to Qv — 0.64. It is interesting to note that at R ~ 1700km 
corresponding to the radius of the Moon, we still obtain rather larger value, Qv = 0.45. 

In Figure H~5l the overlap functions for three representative cases are plotted: R = 600, 1700 and 6400km. As we 
commented earlier, the overlap functions for any radius can be rescaled and become identical if we plot the functions 
against the rescaled variable, y cx fR. Further, for the configuration examined in Figure Q3J the relations 7/(0) = — 1 
and 7y(0) = strictly hold. In those situations, the shape of the noise spectrum shown in Figure [TT1 is the key to 
determine the best value for Qv and the overlap function jv for R ~ 600km eventually becomes the best shape to 
achieve the maximum value of Qv = 0.64. Also, it turned out that the sensitivity for the Moon becomes about three 
times larger than that for the Earth, Qv = 0.15 (see appendix |Pl for comments on detectors on the Moon). 
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V. SEPARATION 



As discussed so far, a simple analysis with correlation signal C a b of two detectors allows us to detect a mixture of I- 
and T^-mode signals, but we cannot extract each of them separately. In order to disentangle these two signals, in this 
section, we discuss the problem of the I- and T^-mode separation considered in Ref. Q • After describing the simplest 
case using a set of four detectors in Sec IV A[ we generalize it to the multiple-detector case in Sec. IV Bl and present 
a statistical framework to achieve the optimal sensitivity. Based on these theoretical backgrounds, in Sec. IV Ci the 
correlation analysis with network of ground-based detectors are examined and the optimal values of SNR are derived 
for each I and V mode. 



A. Analysis with two correlation signals 



Let us begin by considering the simplest case that two pairs of interferometers (a, b) and (c, d) are available. We 
write down their correlation signals as 



Cab(f) 



8tt 



{ll,ab{f)l{f) + 7V,afc(/)V(/)}, C cd (f) 



-{nMfW) + 7vMf)v(f)}. 



(45) 



Making their linear combination, / or V mode can be removed, and a pure V or I mode is separately extracted. 
Except for trivial scaling, the unique choices are 



lV,abC c d — IV.cdCab 



&7T &7T 

I{lI,cdlV,ab ~ JI,ab7V,cd), JI,abC c d ~ H ,cdP a b = ~=-V (jv,cdl I ,ab ~ lV,abl I ,cd) ■ (46) 



Meanwhile, the noise spectra for these combinations are proportional to 

JV,abNcd + Jv,cdNab, 7j ,ab^cd + ll,cd^ab, 



(47) 



where quantity N a b indicates the product of noise spectra, M a b = N a (f)Nb(f). Taking account of proportional factors, 
the broadband signal-to-noise ratio for a pure I or V mode becomes 



SNRf = 



SNRf 



/16tt 
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16tt 
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I 2 {lV,cdll,ab ~ lV,abll,cd) 2 
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V 2 (jV,cd1l,ab ~ lV,abll,cd) 2 



For detectors with identical noise spectra with Af(f) = 7V(/) 2 , the compiled overlap functions are defined as 



I,ab:cd 



H,cdlV,ab — lUablV.cd 

+ nfl U72 : 

UV,ab ^ IV,cd\ 



r 



V.ab:cd 



_ 7 V 2 c<ffL 1 ab — lV,ablI,cd 
U2 , 2 11/2 
Ul,ab ^ II,cd\ 



(48) 
(49) 

(50) 



With these functions, the broadband SNRs for the separated two modes can be estimated from equation (|39]l just 
replacing the term 7// + jvV with Tj^-dl and Tv,ab-.cdV. In this sense, the complied overlap functions represent 
the sensitivities to the / and V modes after the separation. 

In Figure fT6|. as a specific example for the mode separation, we consider the HL and CL pairs and plot the compiled 
overlap function, as well as the overlap functions for each pair. Although the shapes of the functions T a b :cc i are very 
complicated, resultant values of the normalized SNR estimated from equations (|4"5|) and (|49[) are 0.11 for the V mode 
and 0.31 for the / mode, which are very close to the values presented in Sec. IIV Bl (Qv = 0.11 for CL, Qj = 0.32 for 
HL). Note also that for other combinations, the / and I^-mode separation can be performed with nominal changes to 
the naively expected sensitivities Qji t v},ab- This will be discussed in details in Sec. IV Cl 



B. Analysis with multiple data set 

Next consider the generalization of the previous analysis to the cases with multiple data set. For a network of n 
detectors, we can make totally n t = n(n — l)/2 independent correlation signals 

Ci = ^(' m I + 'rviV) (51) 
5 
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FIG. 16: Overlap functions for the unpolarized / mode (dashed curves), and the circularly polarized V mode (solid curves). 
The upper panel shows the results for the Hanford-Livingston (HL) pair. The middle one is results for the LCGT-Livingston 
(CL) pair. The normalized SNRs Si t v (with the adv LIGO noise spectrum) are also presented. The bottom one show the 
compiled functions r^y (Eq.([50}) made from both pairs. 



with i = 1, ..,n t . Here, we use the single suffix i to represent a detector pair for which we have assigned two suffixes 
so far, such as ab in equation (|45p . For the number of detectors with n > 2, the number of output signals becomes 
n t > 3, and this implies that we must deal with the over-determined problem in order to separate a mixture of /- and 
y-modes, because the number of observables exceeds the number of target parameters, / and V. In what follows, we 
will discuss the signal-to-noise ratios expected from the optimal data analysis. 

Let us examine a straightforward extension of the analysis in previous subsection. Provided the original data set 
of correlation signals, {Cj}, obtained from all possible pairs of detectors, we can make the linear combinations, Dn 
and Dyi, which respectively eliminate the variable V and /: 

D n = dul, D Vi = dviV, (52) 

where dn and dyi denote some numerical coefficients, appropriately chosen for removing the contribution from V and 
/, respectively. Note that the number of independent combinations labeled as i is n t — 1 3 . The associated covariance 
matrices for intrinsic noises, -Mj^ and My-, are expressed in terms of the quantities, jn, ■fvi and A/i (noise spectra). 
Then, the resultant total SNRs for the optimal combinations of data sets {Dn, Dyi} are 

SNRj cx d Il M~ 1 d I]1 SNRy cx d Vi M. v \A V j. (53) 



3 One simple example is to make Dn = d - (yvi/jVnt) C« t = {yn - 1in t lVihvn t ) I and D Vi = d - (yn /~/i„ t ) C nt = (yvi - 
IVntlli/jIn^V. 
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For rit 



3, the above expressions can be recast in a rather simple form: 



SNR/ oc I 



SNRy oc V 




(54) 



(55) 



Equations (|54|) and (|55|) are symmetric with respect to the suffix i and they do not depend on the specific choice of 
the data sets {Da, Dyi}- Self-consistently, the above expressions recover the previous results, (|48|) and (|49|) . if we 
set 7/3 and 7^3 to zero. Indeed, the symmetric expressions ([54| and (|55|) for n t = 3 generally hold for the cases with 
rit > 3 and we will use these forms to estimate the SNRs for optimal combination of five ground-based detectors. In 
appendix[Cl a brief sketch to derive the symmetric expressions for the nt > 3 cases is presented. Multiplying the factor 
2(16ir /5) 2 df and integrating over entire frequency range, the narrow band SNRs (|54l) and (|55|) can be generalized to 
the broadband SNRs 
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(56) 



(57) 



Similar to the one defined in previous subsection, we define the effective overlap functions for detectors with identical 
noise spectra, which represent the optimal sensitivities to the / and V modes: 
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For sensitivity only for the / or V mode, we also define 



off, 70 



1/2 



(58) 



(59) 



which correspond to the effective overlap function (with identical noise spectrum) for the traditional analysis in the 
absence of V or / mode. In the following, we use the notation Qjq for the normalized SNR with effective function 
r e ffjo- Based on these definitions, the ratio R is given by 



R = 



i 7i/7 



iv 



off./o 1 cffyo 



(60) 



As increasing the number of detectors, the functions Teff^o and ^ e s,vo monotonically increase. For a large numbers 
of detectors, however, the numerator y\ ' jn"fiv can be regarded as a summation of random numbers, and the ratio 
R is expected to decrease quickly. We will see this numerically in next subsection. 



C. Optimal SNRs from ground-based network 



We are in position to evaluate the broad band SNRs for optimal combination of network of five detectors, A,C,H,L 
and V. For networks made by three detectors among five detectors, there are 5C3 = 10 possible networks. In the 
same way, we can make 5C4 = 5 networks for combinations of four detectors. Numerical results for detector networks 
are presented in Table HVl In Figure [171 we also provide the normalized SNRs for various combinations of detectors, 
showing the overall behaviors against the number of detectors. 

^,From Table HVl and Figure [T7l several diagnostic features are summarized: 

(i) To realize a good sensitivity to the / mode, the HL pair has a crucial role. This is due to their small separation. 
Without the pair, we have at most SNR/ = 0.20 and SNR/o = 0.21. Including the two detectors, the value 
SNR/ becomes more than 0.3. 
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TABLE IV: Normalized SNRs Qi, Qy and Qio for network of detectors. 
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FIG. 17: The normalized SNRs Qi (left panel; filled triangles), Qio (left panel ; open circles), and Qv (right panel; filled 
triangle). The horizontal axis is the number of detectors. We slightly shift the points for Qio to the right. 



(ii) The combination AHL does not have a good sensitivity to the V mode and we obtain Qy = 0.04. This is because 

the orientation of AIGO detector is specialized to achieve the best sensitivity to / mode in combination with 
LIGO detectors. In fact, they are aligned to have large overlaps (in relation to (i)). Nevertheless, the sensitivity 
to the V mode can be improved by adding the LCGT or Virgo detector. In contrast, the value Qio increases 
only 10% even if we increase the network from AHL to ACHLV. 

(iii) Comparing SNR/ with SNR/o, we deduce a tiny amount of statistical loss for the sensitivity to the I mode, 
caused by adding a new target parameter, V. It is not preferable to get a significant statistical loss by dealing 
with the circular polarization mode whose fraction is naively expected to be small. But Table IIVI indicates 
that even when the V mode is added as observational targets, detection efficiency for / mode remains almost 



21 




50 100 150 200 50 100 150 200 

f [Hz] f [Hz] 

FIG. 18: Left: effective overlap functions, F 2 g 7 (solid), F 2 S ^(dotted), and T^g J0 (long-dashed). Right: same as in left panel, 
but for the network of five detectors. Two curves for T 2 B I and r^ ff I0 are nearly overlapped. 
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FIG. 19: The function R 2 for the CHV (dotted curve) and the ACHLV networks (solid curve). 



unchanged. For three-detector networks, the relative loss is largest for CHV and reach about 24%. Increasing 
the number of detectors, the maximum loss is reduced to 4% for four-detector network, and further reduced to 
1% for five-detector network. 

To further give a deep insight into the diagnosis (iii), in Figure fTSl we plot the effective overlap functions r c ffj, T e sjo 
and T c sy for specific networks of CHV (left) and ACHLV (right). For network of CHV, there exists characteristic 
pattern at / > 40Hz with period A/ = 18Hz. The main reason of this comes from the fact that all the pairs, CH, 
CV and HV, have the separation angle (3 ~ 80°, leading to the frequency A/ = c/(ARe sin(80°/2)) ~ 18Hz. Focusing 
on the differences between T c f[j and T c gjQ, we find that while the differences are manifest at low-frequency in CHV 
system, the functions Teg,/ and r e ffjo in ACHLV system become almost identical even at low- frequency. This is 
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clearly quantified if we plot the ratio R defined in equation (|60|) . Figure [19] reveals that the magnitude of the ratio R 
is significantly reduced for the ACHLV system, suggesting the fact that the statistical loss is negligibly small. In this 
respect, negligible statistical loss may be another merit for a network with a large number of detectors. 

Finally, from numerical results given above and with a help of equation (|42[) . we summarize the signal-to- noise 
ratios, SNR/ and SNRq (not normalized ones) with noise spectrum of advanced LIGO. For the five-detector network, 
assuming the flat spectra Slew = const and SIqwII = const, we have 

SNR, = 1.64 (O^o) (^) ~ 1/2 , SNRy = 0.749 (^^) (^) . (61) 



VI. SUMMARY 



In this paper we present prospects for measuring the Stokes V parameter of stochastic gravitational wave back- 
grounds via the correlation analysis. This parameter characterizes the asymmetry of amplitudes of right- and left- 
handed waves. As the parity transformation interchanges the two polarization modes, it can be regarded as the basic 
observational measure to probe parity violation. We made detailed analyses for the basic properties of the overlap 
functions 7/ and 7y, especially their dependencies on geometry of detector configurations. 

In contrast to studies only for the unpolarized / mode (equivalently, the energy spectrum Oqw), we need to develop 
a new statistical framework to deal with rich structures caused by multi-dimensionality of target parameters. We 
provide an optimal method that will be applicable to various problems of gravitational- wave backgrounds. Based on 
our new method, we estimated sensitivities of the planned and proposed next-generation interferometers to the V 
modes. We found that it is important to have a large number of detectors in order to reduce possible effects due to 
correlation between target parameters. 

We would like to thank M. Ando, N. Kanda and M. Ohashi for useful conversations. This work was in part supported 
by a Grant-in- Aid for Scientific Research from the Japan Society for the Promotion of Science (No. 18740132). 



APPENDIX A: TENSORIAL EXPANSION 



In this appendix we present tensorial decompositions of the overlap functions 7/ y ( see a l so Ref.[Hj]) defined for 
pair of detectors a and b at positions x a and xi>. They are expressed as 

7/,o*(/) = ^ dn [{F+F+* + e^ nm ] , (Al) 

and 

IvM = f dn[i {F+FF - F*Fj}~*} e iynrn ] , (A2) 

with x a — Xf, = Dm (D: distance, m:unit vector) and y = 2irfD/c. The beam-pattern functions F^ are written by 
the the polarization tensor e p and the detector tensor d a as 

F^d a :e p (n)=d^. (A3) 

Here the detector tensor d a is given by two orthonormal vectors u a and v a as d a = (u a ®u a — v a ®v a )/2. Therefore, 
the overlap functions are formally written as 

n(f) = r /ltfW (/)44i. ivU) = Tv^MmAi- ( A4 ) 

In these expressions we defined 

IW/) = 2. J dn [e+(n)e+(n) + e*.(n)e fe »] e^ nm (A5) 

and 

ry,««(/) = ~ J dn [e+(n)e* (n) - e*(n)e+ (n)] e^ nm . (A6) 
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There are apparent symmetries with respect to the subscripts of the tensor Ti^jki- For example, it is invariant 
under the replacement i *-+ j or «-> (k,l). Furthermore, with using the correspondences e^j(— n) = e^t-(n) and 

e^(— n) = — e^-(n) for parity transformation, the tensor Yi^jki is a real function taking a same value at m and — m. 
^From these symmetries, the tensor Tx t ijki{f) is given by a combination of basic tensors ra, L and &y as follows; 

Fi,ijkl{f) = inSijSki + a I2 (Si k Sji + SuSjk) + a IZ {8 tj m k mi + Skim z mj) + a IA m i m : jm k mi 

+a I5 (8 ik m : jmi + 5 jk m t mi + S 3 imim k + 5am jim*,), (A7) 

with the expansion coefficients an. These coefficients are given in the following manner. We firstly fix the direction 
vector m = (f , 0, 0) and calculate the components Ti t ijki(f) f° r ^, = ( x j x i x i x )i ( x i x i V> 2/)> ( x > Vi x i v\ (Vi Ui V) 
and (y, y, z, z). Then we can solve the coefficients an with their five independent combinations. After some calculation, 
we obtain them in terms of spherical Bessel functions with argument y as 

fO f 20 5 5 15 5 , An , 

an = 0, a I2 = j - —j 2 + a/3 = — j32 - a u = -j 4 , a I5 = —j 2 - — j 4 - (A8) 

A detector tensor d a is usually traceless {dfjSij = 0), as we measure quadrupole deformation of space e.g with 
interfering laser beams of two arms. Then the first and third terms in equation (|A7|) do not provide contribution to 
the overlap function 7/. With angular parameters (/3, <7i,cr 2 ) for a given detector pair (a, b) on a sphere with radius 
R (see figure 1), we can set their positions as x a = i?(cos/3/2, 0, sin/3/2) and x b = i?(cos/3/2, 0, — sin/3/2), since only 
their relative positions are relevant. In this case we have m = (0,0,-1), and the two unit vectors u a and v a are 
written by 

u a = cos(7i(sin/3/2, 0, — cos/3/2) + sin<7i(0, 1,0), v a = — sincri(sin/3/2, 0, — cos/3/2) + coseri(0, 1,0), (A9) 
while we can put 

u b = coscr 2 (sin/3/2,0, - cos/3/2) + sincr 2 (0, 1,0), v b = -sincr 2 (- sin/3/2, 0, - cos/3/2) + coser 2 (0, 1,0) (A10) 

for the second detector b. With plugging in these expressions into equation (|A4|) we obtain equation (|22"|) . 

Similarly, the tensor Tvijkt is invariant with replacements such as i <-> j, but it is asymmetric for the replacement 
<-> (k, I) or Til 2 — > — TOj. The tensor is real due to the parity relation as for Tjjj k i. We found that it is expanded 
with the basic tensors as 

^Yjki = avi(^ikS 3 i + LOtiS-jt, + Wjk6u + LOjiS ik ) + ay 2 (ui k 7n 7rn + wumjTrik + Lo jk mimi + uj 3 imim k ) (All) 
with ujij = eij k m k (eij k : antisymmetric tensor). In this case, the expansion coefficients ayi are solved as 

1 5 

avi = ji ~ av2 = tJ3- (A12) 

We can derive equation (j2l)|) as in the case for the / mode analyzed above. 

When two detectors a and b are on a same plane, the tensor cjy cannot have component with respect to the two 
dimensional projected space to the plane, and we have identically jv = with equations (|A4p and (|A11[) . 

Note that the tensor Tj is an even function of rrii but the tensor Ty is an odd function reflecting its handedness. The 
function 7/ is given by spherical Bessel functions ji with even i, while the function 7y is with odd i. The asymptotic 
behaviors of the spherical Bessel functions are 

/ n 1 ( (n+ l)ir\ , , , 

Uv) ~ - cos [y - 2 J (A13) 

at y — > 00. In the same manner the peaks of the function |7/| are at y ~ (N + 1/2)tv (N: natural number) and those 
for \"fv\ are at y ~ Ntt. Therefore the zero points for jv a- n d 77 are offset by Ay = tt/2 at large y. 
At the low frequency limit y — > we have 

MV) - T^yyy, (A14) 

and the asymptotic behaviors of the overlap functions are 7y — > and 7/ — > dydy/2. For two traceless tensors rib- 
and d^, the combination d1jd\j is the unique scalar quantity written by their tensor product. In terms of the angular 
parameters in the main text, this limit is given by cos 4 (/3/2) cos(4A). 
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APPENDIX B: PROBABILITY DISTRIBUTION FUNCTIONS FOR CORRELATION ANALYSIS 

In this appendix, we study probability distribution functions associated with correlation analysis, following Ref . (27j . 
With Fourier space representation, each data stream s a (f) is made by gravitational wave signal H a (f) and noise n a (f) 
as 



s a (f) = H a (f)+n a (f), 



(Bl) 



and we define its noise spectrum 



K(f)n a {f')) = -S D (f-f)N a (f). 



(B2) 



We assume that the correlation between noises are negligible (namely {n* a n b ) = for a ^ i), and the amplitude of the 
signal (H a (f)*Hb(f)) is much smaller than that of the noise (n a (f)*n a (f)). These are the conditions where correlation 
analysis becomes very powerful. We divide the positive Fourier space into frequency segments F v (i = 1, ...,N) with 
their center frequencies {f v } and widths {5f v }. In each segment the width 5f v is much smaller than /„, and the relevant 
quantities (e.g. Ogw(/)j Jab(f)) are almost constant. But the width is much larger than the frequency resolution 
Af = T~^ s (T D b s : observation period) so that each segment contains Fourier modes as many as Sf v /Af » 1. 

For correlation analysis we compress the observational data sj(/) by summing up the products s*(f)s b (f) (a ^ b) 
in each segment F v as 



fj, v = ^ S *a(f) S b(f), 

feF v 



(B3) 



where we omitted the apparent subscript {ab} for the compressed data {fJ. v } for notational simplicity. As the noises are 
assumed to be uncorrelated, the statistical mean (fi v ) is caused by gravitational wave signal. After some calculations, 
we have a real value 



(aO = £ (H a (fyH b (f)) 



yU7aM + V-f ab y) — . 



(B4) 



The fluctuations around the mean (fj, v ) are dominated by the noise under our weak signal approximation, and its 
variance a% for the real part of ju, v becomes 



a 2 v = N a (fv)N b (f v ) 



Sfv 

8Af 



(B5) 



As the number of Fourier modes df v /Af in each segment is much larger than unity, the probability distribution 
function (PDF) for the real part of the measured value \i v is close to Gaussian distribution due to the central-limit 
theorem as 



p(Re[fi v }) 



1 



y/27 



exp 



(Rc[/i„] - (/i^)) 2 
2ul 



(B6) 



Here we neglected the prior information of the spectrum Uqw(/) and n(/). From equations (|B4[) and (|B5[) , the signal 
to noise ratio of each segment becomes 



SNR, 



16tt 
"1T 



obs 



26f v 



(IjabJ + Vjaby) 2 



Na(f)N b (f) . 

Summing up the all the segments quadratically, we get the total signal to noise ratio 



(B7) 



SNR' 



16tt 
"1T 



T, 



obs 



df 



(Ijab.I + Vlab,vf 
Na(f)N b (f) 



(B8) 



This expression does not depend on the details of the segmentation {F v }. The same results can be derived by 
introducing the optimal filter for the product N*(f)N b (f) to get the highest signal to noise ratio (see e.g., Ref.p^j). 
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APPENDIX C: DERIVATION OF OPTIMAL SNRS FOR MULTIPLE DETECTORS 



The derivation outlined in Scc lVBl is intuitive, but it is algebraically complicated to derive the final expressions. For 
example, we need to deal with large-dimensional noise matrices M-nj and Mvij that have off-diagonal components. 
In this Appendix, we make a simple explanation for the structure of equations (|54|) and (|55[) . and derive useful 
expressions valid for arbitrary number of detectors n. Based on Appendix [B] we consider the summed correlation 
signals {/i^} in a fixed small band F v with its bandwidth 6 fa as in equation (B3). Later, we will sum up all the bands 
to get the total SNRs. 

In actual observation, we cannot exactly measure the expectation values 



8ir r 

—{~f U I + ~f Vi V}b v . 



(CI) 



Rather, the measured values {n' vi } fluctuate around the true values Cib v with variances N 2 b v . Here the ratio b v — -^j 

is the number of frequency bin in the band /„ with the frequency resolution A/ = TI7 . The multi-dimensional 
probability distribution function P{{n' vi }) has a form 



with the kernel 



P(KJ) « exp [-K] 

(M$n - (MOT)) 2 



a/? 



b v . 



(C2) 



(C3) 



The structure of the distribution function P(I e , V e ) for the estimated values {I e , V e } is obtained by the replacement 

8tt 



Mi 



-(llile + lViV e )b v 



(C4) 



Since the expression in the large parenthesis [• ■ ■ ] in equation (|C2[) becomes a quadratic function of the target 
parameters {I e , V e }, their expectation values are {J, V}, and their covariance noises matrix is proportional to 



an aiv 
a iv a vv 



1 



anavv — a>iv 



avv —aiv 
—aiv a n 



(C5) 



with the matrix elements, an = ^ 4 7|j/M, avv = J2i Ivt/Ni and aiv — Y^ixlulVi) /-M- Note that the off-diagonal 
element aiv has information for the statistical correlation between I- and V-modes. We define the ratio R by 



R 



aiv 



yfanavv 



The ratio of the expectation values squared {I 2 , V 2 } to the variances of their noises are proportional to 



2 anayy — a IV 
avv 



V 



■ aiiayy — Ojv 
an 



(C6) 



(C7) 



These expressions exactly coincide with equations ([51)1 and (|55[) in the case of n t — 3. Indeed, using Mathematica, 
we confirmed that the above expressions faithfully reproduced the same results as obtained from the direct estimation 
with equation (|53p . up to n t — 8. By summing up all the frequency bands and properly dealing with the number of 
bins b v = Sf v /Af (as in Appendix IB"]), we can easily evaluate the broadband SNRs for large number of detectors with 
the following the expressions: 



SNR 2 



SNRy = 



(T) 

16tt 
IT 



-M3bs 



2obs 



dfi 



2 anavv 



l iv 



dfV 



avv 
2 anavv — a 2 v 



an 



For analysis only with the J-mode (or Qgw), the broadband SNR for the I mode is evaluated by 

2 



SNR' 



10 



/16?r 



dfl 2 an 



(C8) 
(C9) 

(CIO) 
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We use this expression as a reference to analyze effects caused by estimation of multiple parameters. We can express 
SNR/ as 



SNRj 



16tt 
~5~ 



2~obs 



p CO 

2 / dfl 2 au{l - R 2 
Jo 



(Cll) 



Therefore the ratio R characterizes the loss of SNRs due to the increase of the number of observable parameters. A 
similar argument holds for the V mode. 



APPENDIX D: DETECTORS ON THE MOON 



It has been discussed that, on the surface of the Moon, the high- vacuum level and rich three-dimensional surface 
structure are suitable to build a gravitational wave interferometer with very long armlength (see e.g., Ref.j29j). 
Meanwhile, the north and south poles of the Moon seem to be preferable and are thought to be special places for 
human activities as well as astronomical observation (see e.g., Ref . [30| ) . The rotation axis of the Moon is nearly 
perpendicular to the ecliptic plane. Therefore, at the rims of craters near the two poles, the sun-light is available 
during most of one Moon's day (~30 Earth days). In contrast, the bottom of craters around the poles is at permanent 
night with exceptionally stable temperature environment around 40K, and we might obtain trapped water ice, from 
which we can produce hydrogen and oxygen (fuel for rocket engine) with electrical decomposition. Away from the 
pole areas, the surface of the Moon has severe physical conditions with a large temperature variation typically from 
~ 100K (night) to ~ 400K (daytime). Therefore, when we build detectors on the Moon, location near the poles would 
be the most natural choice. 
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